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Abstract 



Based on the Nakajima-Zubarev type nonequilibrium density operator, we derive a hyperbolic 
hydrodynamical equation. Microscopic Kubo-formulas for all coefficients in the hyperbolic hy- 
drodynamics are obtained. Coefficients a^s and /3j's in the Israel-Stewart equation are given as 
current-weighted correlation lengths which are to be calculated in statistical mechanics. 
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I. INTRODUCTION 



The hydro dynamic model is one of the widely applied phenomenological models for rel- 
ativistic heavy ion collisions. According to the recent detailed analyses, weak but non- 
vanishing viscosity is significant for the qualitative discussion of vi at RHIC Because 
the Navier-Stokes equation is a parabolic type equation, naive relativistic extension is not 
consistent with causality, and as a result, the numerical solution becomes unstable fl. The 
relativistic causal hydrodynamic model which is a hyperbolic type had been introduced 
phenomenologically by Israel and Stewart more than thirty years ago 0, Q]. In the past 
decades, many proposals appear for the appropriate equations of the relativistic causal hy- 
drodynamic model Q-7|. Most of the works seem to belong to the semi-classical phenomeno- 
logical approach based on the Boltzmann equation. In this paper we derive a hyperbolic 
type hydrodynamics as a second order hydrodynamics based on a nonequilibrium density 



operator method [8( . Romatschke and Moore et al. also discussed hydrodynamics based on 
a nonequilibrium entropy density (9), 1(3]. In the present paper, more explicit formulas for 
the coefficients are issued. 

The hydrodynamic model is a phenomenological model of the macroscopic point of view 
which is based on the coarse graining. The total system is constructed as a patchwork 
of a local system corresponding the statistical systems in equilibrium. The local system 
at each space-time point is expected to be microscopically large but macroscopically one 
point and moving with four velocity U^ 1 . Thermodynamical quantities and the coefficients 
in the hydrodynamical equation are the local quantities as the functions of position through 
the local thermodynamical parameters such as temperature T(x) and chemical potential 
n(x). The macroscopic model based on the coarse graining is justified by the fact that the 
microscopic correlation length is much smaller than the typical scale of the macroscopic 
dynamics. More than fifty years ago, Iso et al. discussed the applicability condition of the 
hydrodynamic model for the multiple production in high energy pp collisions based on the 
comparison between the length of the canonical correlation of the current operators and the 



scale of the change in the solution of the hydrodynamical equation [11 ]. 

Regarding the self-consistent macroscopic model, before we start to solve the hydrody- 
namical equation, all functional forms of the coefficients should be fixed. All coefficients of 
the hydrodynamical equation are given as functions of thermodynamical parameters, T and 



/i. The derivation of them is a major task for statistical mechanics. The well-known Kubo 
formulas for the viscosities and heat conductivity provide us ways to evaluate the coefficients 
of the Navier-Stokes equation by statistical mechanics in equilibrium. However, the method 
to calculate the additional coefficients in the hyperbolic hydrodynamical equation, the so 
called Israel- Stewart equation, have not yet been established. As far as the authors know, 
only the kinetic calculation by using the Grad 14 moments method is adopted. The aim of 
this paper is to establish a way to calculate all coefficients in the hyperbolic hydrodynam- 
ics as the functions of temperature and chemical potential. Based on the nonequilibrium 
density operator method, we systematically derive microscopic formulas for the coefficients 
ccj's and /3j's in the Israel- Stewart equation. All coefficients are expressed as the current- 
weighted length of canonical correlation which are to be calculated in statistical mechanics 
in equilibrium with T and fj,. In general, evaluation of the canonical correlation is not 
easy, but our formulas can be good targets for the future Lattice QCD simulations and the 
hadro-molecular calculation. 



II. HYDRODYNAMIC MODEL 

Hydrodynamical equation is composed of the energy-momentum conservation law 

<9%, = (1) 

and the charge conservation law 

d% = 0. (2) 

The flow of the fluid is described by a four velocity which is a normalized time-like vector, 
U^U^ = 1. By using f/ M , the space-like projection operator is defined as = g^ v — U^U U . 
With £/ M and A MI/ , we can define four projection operators of the second rank symmetric 
tensor: 

(4) 

^HPjjujja + U^ljc^vp + ^JJVJJP} ; (5) 

(6) 
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T( At!y | pfJ ) stands for the traceless part of the spatial component, S{^\ pa ) stands for the trace 
of the spatial component, U( tll '\ pcr ) stands for the time-like time-like component, and V( pu \ pcr ) 
stands for the time-like space- like or the space-like time- like component, respectively. These 
operators form a complete set 

s{ pu \ pa ) + rc^H + Ki^D + VH pff ) = - {g w g va + g^g up ) = i. (7) 

jjh corresponds to the direction of the time on the local rest frame at x p and A Mt/ specifies 
the space-like coordinate, respectively. Hence, on the local rest frame, becomes (1, 0, 0, 0) 
and becomes —8 l \ respectively. Throughout this paper, Greek indexes, fi,u,p,a etc. 
stand for the Lorentzian (0,1,2,3), and Latin indexes k etc. stand for the Euclidian (1,2,3), 
respectively. A space-time point (x, t) is denoted as x M for the abbreviation, r = x^U^ stands 
for the time on the local rest frame (comoving frame of the fluid element), and D = U^d^ 
stands for the derivative by r, respectively. 

Energy density is given as the time-like time-like component of the energy-momentum 
tensor and pressure is given as the the trace of the spatial part, respectively, 

e = T^U^Uu, (8) 
P = ~T^A^. (9) 

Once Lagrangian of the system is given, the energy-momentum tensor T pu is defined as 
the space-time translation operator. However, the definition of the flow U p of a relativistic 
fluid is not trivial. Eckart (l2fl and Namiki et al. [13j define based on the particle flow 

n 

(P-frame) and Weinberg used charge current [14j . In the textbook by Landau and Lifshitz, 



U p is defined as a eigenvector of T^ u 15j. Recently another possible choice is proposed by 



Tsumura et al. 



16| . In this paper we adopt the Landau-Lifshitz frame (E-frame) where 



is an eigenvector of the energy momentum tensor: 

T» U U U . = eU p . (10) 
Stress-shear tensor is given as the spatial traceless part of the energy momentum tensor: 

n pu = m p(7 )T pa . (11) 
Using the projection operators, we can denote energy density as: 

e = U^U v UrU)T pa , (12) 



and pressure as: 

P = l -/\, v Sr\ P «)T p % (13) 

respectively. The time-like space-like components of the energy momentum tensor should 
correspond to the heat flow, however, because of Eq. (fTOj) these components must vanish: 

v( M T CT )^ pCT = 0, (14) 

on the E-frame. 

The charge density n and the charge current J M are given by n = J^U^ as usual, however, 
is proportional to the energy flow on the E-frame, and the charge current contains J M 
which is perpendicular to U^: 

J» = nW + J**, (15) 
in corresponds to the heat flow on the E-frame [15]. 



III. NONEQUILIBRIUM DENSITY OPERATOR 



The nonequilibrium density operator is given as 17H19] 



p = Q- L exp(-A + B) 



with A being: 



B being: 



+Ux' ) t')& v (f3(x' 1 t')^x' ) t')) 



U^x-x')^ > 



(17) 



A = JdS» {/3(x, t)U»(x, t)% v {x, t) - P(x, t)fi(x, t)j„(x, t) } , (18) 



(19) 



and with Q being normalization Q = tr[exp(— A + B)]. The expectation value of an operator 
0(x) is given by (0(x)) = tr pO(x) . 

A corresponds to the local equilibrium density operator, /)/ = Q^exp^—A), with Qj l = 
tr exp(— A) being normalization. dS^ is a three-dimensional hypersurface which locally 
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corresponds to the comoving volume element of the fluid. The so-called matching conditions, 
(e) = (e)i and (h) = (n)i, are imposed for /3(x,t) and /3(x,t)fjL(x,t) in pi. If the expectation 
value of some operator O(x^): 

(O(xn)i = QThrfadia?)], (20) 

is a local quantity in the macroscopic sense, on the local rest frame at that point, the 
expectation value is the same to the one in the corresponding equilibrium state where the 
inverse temperature is 1/T = 0(x,t) and the chemical potential is p(x,t), respectively. 

B stands for the influence of the thermodynamical forces such as the disturbances in the 
flow and the thermodynamical parameters. The adiabatic limit ( — > 0+ should be taken at 
the final stage of the calculation. 

According to the standard procedure of the linear response theory, we expand the nonequi- 
librium density operator p at around the local equilibrium operator p\ up to the linear term 
of the thermodynamical forces. Expectation value (0(x)) is given as a canonical correlation 
with the force term: 

«*,)> = <0 W>1 + .„„ /,Ve-<-». (o W ,^<) W )^<))) 

+ Mm J dVe-^*-*'^ (6(x), J p {x')d' p {(3{x')p{x'))) , (21) 

with (Oi,0 2 ) being: 

{Ox, 6 2 ) = (6i f 1 d\e XA d 2 e- XA ) l - (di)t(d 2 )i. (22) 
J o 

IV. DERIVATIVE EXPANSION WITH THE THERMODYNAMICAL FORCES 

If the change of the thermodynamical parameters and the flow are slow enough and can 
be treated as almost constant during the microscopic relaxation length which is determined 
by the canonical correlation Eq. (T2"2"j) . we may adopt the Taylor expansion of the thermody- 
namical forces d' p (f3(x')U a (x ; )) and d' p (P(x')p(x')) at around the position x of the operator 
0(x), and put the forces outside of the integration in Eq. (T2TT) . Then we can obtain expec- 
tation values of T^ v and J M in the power series of the gradient of the thermodynamical 
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parameters and the flow: 

(f^(x)) = (f^x^+fj d 3 x' y* dt'e-^if^ix),^'))^ d p (f3U°) 
+ d 3 x' J dt'e-^ t - t '\f^{x)J p {x')) ) jd p {^) 

+ ( I d 3 x' [ dt'e-^if^ix), (x x ' - x x )J p (x'))) d x d» {fin) , (23) 



(>(x)> = (>(ar)),+ ( y dVy dt'e-^^'H^Cx), J p (x'))J ^ 

+ ( [ d 3 x' J* dt'e-^ t - t '\r(x),(x x ' -x x )J p (x'))) d x d p ((5 l i) 

d 3 x' f Me-«t-V(P(x), (x x ' - x x )f pa (x'))) d x d p (f3U°) . (24) 



Here all coefficients are given as the integral of the canonical correlation with the local equi- 
librium density operator f !22p . If the non- vanishing region of the canonical correlation in 
the integrand are limited only in the microscopic scale, we can regard them as local quan- 
tities in the macroscopic sense and we may substitute a statistical mechanical calculation 
in equiliblium with the corresponding T and \i for the expectation with p\ in Eqs. ( 123]) and 

(El. 

The first terms in the right hand side of (123]) and (124"]) . which are free from the thermody- 
namical forces, correspond to the perfect fluid part: (T^x)); = ((e) i + (P)i)U^U u — (P)^" 
and (J^(x))i = (n)z£/ M . The second terms in the right hand side of (123]) and (124"]) . which are 
the the zeroth order gradient terms of thermodynamical forces, correspond to the Navier- 
Stokes equation. The Kubo-formulas for viscosities and heat conductivity are obtained as 
the coefficients in these terms 



13, 



18] . The third terms in the right hand side of ( 123]) and 



24]) . which are the the first order gradient terms of thermodynamical forces, correspond to 



the hyperbolic hydrodynamical equation 

Decomposing into the space-like components and the time-like components, we can 



rewrite B as: 



B = Km J d A x'e- cu ^ x - x '^ {/3(a/, t')^ v {x! t')&»{U v {x\t')) 
-i v {x',t')d" / (/3(x',t , )fx{x , ,t')) 

-/3(x',t')P'{x;t')d^{x',t')}, (25) 



where: 



By virtue of the matching condition and the E-frame condition (ITUl) . S{f xv \ pa )T ( " J and 
V(^\ pa )TP a don't appear in B. 

Let us denote integration of the canonical correlation in Eqs. (I23p . 0241) as: 

«Oi|0 2 )) = Kin y dV^' rft'e-^-^ (d 1 ( S B,*),0 2 (a:',O) , (27) 

for the abbreviation. Assuming isotropy of the local rest system, we can apply Curie's 
theorem for ((Ojj02)). ((O1IO2)) depends on (x,t) only through the temperature and the 
chemical potential. Because of the translational invariance of the corresponding equilibrium 
system, we can rewrite: 



Hm J d 3 x' J dtfe-tW (di(aj, t), (aP - x»)6 2 {x', t'f) , (28) 
hm + jd 3 xj dte +ct (d 1 (0,0),x^O 2 (x,tf) , (29) 



C- 

{{O x \x»\0 2 )). (30) 



Substituting the projection operator g^ v = U^U U + into equations (1231 and (124|) . and 
applying the Curie's theorem, we can obtain simple forms of the expectation values of the 
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thermodynamical currents: 



(rrn = (tO, + f3((n^n)d p U a - {{^\P))d p {^) 

+P((n^\x x \nn}dxd p U a - ((^\x x \P}}d x d p (^) 

+(3{{^\r\^))Dd p U a - ((nnxx\I p ))A XK d K d p ((3fi), (31) 

(/>) = (/">, - wmwn) - {{p\x p \n)dM^) 

+(3((P\xV T ))d u d p U a - f3{{P\x v \P))d v d p U p 
= (J"), - /3((F|r))A/9 p (^) - «/"|r|/"»I>^(^) 

+/3((/^>^>>A/^ p f/ CT - P((P\x v \P'))K a d c d p U p , (32) 
(P) = (P) l -P((P\P))d„U 1 ' 

-P((P\x»\P))d p d„U v - {{P\x p \P))d p d v {^) 
= (P) l -f3((P\P))dM u 

-(3((P\r\P))Dd v U» - ((P\x p \P))\°d a d v (t3n). (33) 



For the isotropic system, we can factorize scalar coefficients from tensor structures as 
follows: 



((P\P)) = (I\I)5^, (34) 

(<7r^|7r fei >> = (s ik 5 jk + 5 il 5 jk - ^5 kl ^j , (35) 

((n i 3\T\ir kl )) = (W* + 5 il 5 jk - ^5 ij 5 kl ^j , (36) 

((P\r\P)) = (I\t\I)S», (37) 

((P\P)) = (P\P), (38) 

((P\t\P)) = (P\t\P>), (39) 

({P\x j \P)) = (I\x\P)5 lj , (40) 

({7T ij \x k \I 1 )) = L lk 5 jk + S il S jk - \d ij 5 kl \ , (41) 



9 



where scalar coefficients in the right hand side are defined by: 



(n 7r) = 




(42) 


(lT t 7r) = 




(43) 


= 




(44) 


(I\t\I) = 




(45) 


(7T|X|/) = 




(46) 


(J|x|P') = 


^-A^((I,\x u \P')), 


(47) 


(P|x|J) = 


^-A^((P\x u \I,)). 


(48) 



Both in equations fl3T|) and (1321) . because of the isotropy, the first term in the right hand 
side should vanish: 

(tO = 0, (49) 

(ni = o. (50) 

(P)i is static pressure. The second terms in the right hand side of Eqs. (l3~lj) and (1331) stand 
for the Kubo- formulas for the shear viscosity rj s and the bulk viscosity rj v : 



r ]v = f3(P\P'). 



(51) 
(52) 



The coefficient in the second term in the right hand side of Eq. ( 1321) corresponds to the heat 
conductivity k on the E- frame 151 ]: 



, + (P)i , , nr 



(53) 
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Finally, we can rewrite equations (|3T]) .( l32|) and (|33|) as: 

i^n = vsTrr)d P u a + ^-n\t\-n)Trr)Dd p U a 

-(7r\x\I}Trr)d P d„(Pri, (54) 
(h = k f ... (w) f 6 . ) ^ u d v {^) + (I\t\I)A^Dd u (^) 

+f3(I\x\n)Trr)dud p U a , +f3(I\x\P')A^d u d a U r7 , (55) 

(P) - (P) t = -r, v d^ - (3(P\t\P')Dd,U» + (P\x\I)A^d,d u ((3fi). (56) 
Substituting the first order equation (the Navier-Stokes equation): 

\{e)i + {P)iJ 
(P) - (P) l = -r, v d,U^ 

for the second order derivative terms of £7 M and (3fi, we can obtain the Israel-Stewart equa- 
tion: 

Vs 

. 



(i)l+(P)i 

WJ^lrrndA^) + m * in Ai»d v Up) - (P)t) , (58) 

Vs -T]v V ' ' 7 



(P) - (P) t = - Vv d,U» + &^1 D Up) - (P),) + (P|X|7) A^(t). (59) 
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Comparing the above equations ( T57|) . ( T58|) and (|59|) with (2.38) in ref |4|, we can identify 
the aj's and /3j's in the Israel-Stewart equation in our notation. Besides thermodynamical 

11 



quantities and the numerical factors coming from the differences in the definition, in 
ref. [4 1 are current-weighted correlation times; (3q ~ (P|t|P'}, /?i ~ (P|t|P/ and /3 2 ~ (7r|t|7r). 
«j in ref. (4] is cross-current- weighted correlation length; a\ ~ (J|x|7r). As concerns bulk 
current, though the original current operator is P the thermodynamical current in B is P', 
hence, On ~ (P|x|J) in bulk current (2.38a) and — «o ~ (J|x|P') in charge current (2.38b) 
in ref. [4]. 

Normalizing the "state" | tt) , |J) and |P') by using the relation, f l5~Tj) . f }52l and (1531 . we 
define the current-weighted correlation times: 

(/|t|/) 



Tl 



T, 



(n 


t 




(tt| 





[P\t\P') 

and the cross-current weighted correlation distances: 



x„/ 



(60) 
(61) 
(62) 




(63) 
(64) 
(65) 



All these quantities are the functions of the temperature T and the chemical potential /x, 
which are able to be calculated based on the statistical mechanics in equilibrium. 



V. CONCLUDING REMARKS 



Based on the nonequilibrium density operator, we have derived the hyperbolic hydro- 
dynamical equation. All coefficients in the equations are expressed as the integration of 
the canonical correlations of the current operators which are to be calculated in statisti- 
cal mechanics. Our discussion is nothing special for the relativity but we simply expand 
the thermodynamical forces in their derivatives. The additional coefficients in the Israel- 
Stewart equation a^s and /3j's are given as the current-weighted correlation times and the 
cross-current-weighted correlation distances. 
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The current- weighted correlation time corresponds to the relaxation time of the current. 
The evaluation of the relaxation time (3 j is essentially the same to the calculation of the Kubo- 



formulas for the transport coefficient 20]. According to the hadro-molecular simulation, 
each current exhibits it's own relaxation even if the basic dynamics is common and once we 
succeed to figure out the behavior of the relaxation of the currents, we can easily evaluate 



both the relaxation time j3i and the transport coefficient 



2 111. This method has also been 



applied to the calculation of the viscosity of a hadron gas|22 



2311 . a,-'s are the correlation 



distances between different currents which have not yet been investigated. 
As Iso et.al. discussed, the comparison between the microscopic corre 



ation length and 

macroscopic scale in hydrodynamics is the touch-stone of the model [llj. The canonical 
correlations exhibit typical microscopic scales of the system. The evaluation of and /3j or 
quantities in Eqs. (62)~(67) based on the Lattice QCD, or hadro-molecular calculation, will 
provide us the key informations to justify the hydrodynamic model of the QCD matter. We 
have treated only the linear term of the thermodymical forces, but higher order purturbation 
or non-perturbative treatment would be appreciated in the microscopic calculation of the 
canonical correlation. 

The bulk current, Eq. ( J26i) is a natural extension of the P' in jl9| and 18 1 where the 



charge conservation is neglected. If a conserved current exists, P cx e is not enough to make 
P' vanish. Therefore, even if (P) = is achieved in a high energy region, the vanishing 
of the bulk viscosity may not be trivial. 
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